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We use the linear 5 expansion, or optimized perturbation theory, to evaluate the effective potential for the two 
dimensional Gross-Neveu model at finite temperature and density obtaining analytical equations for the critical 
temperature, chemical potential and fermionic mass which include finite A' corrections. Our results seem to 
iinprove over the traditional large-A' predictions. 
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I. INTRODUCTION 



The development of reliable analytical non-perturbative 
techniques to treat problems related to phase transitions in 
quantum chromodynamics (QCD) represents an important do- 
main of research within quantum field theories. The appear- 
ance of large infrared divergences, happening for example in 
massless field theories, like in QCD 1 1], close to critical tem- 
peratures (in field theories displaying a second order phase 
transition or a weakly first order transition can only be 
dealt with in a non-perturbative fashion. Among the analytical 
non-perturbative techniques one of the most used is the 
approximation |3]. Though a powerful resummation method, 
this approximation can quickly become cumbersome after the 
resummation of the first leading contributions, like for the 
N = 3 case which regards QCD. This is due to technical dif- 
ficulties such as the formal resummation of infinite subsets 
of Feynman graphs and their subsequent renormalization. In 
this work we employ an alternative non-perturbative method 
known as the linear 5 expansion (LDE) |4] to investigate the 
breaking and restoration of chiral symmetry within the two 
dimensional Gross-Neveu model at finite temperature (T) 
and chemical potential (/j). As we shall see, the LDE great 
advantage is that the actual selection and evaluation,including 
renormalization, of the relevant contributions are carried out 
in a completely perturbative way. Non-perturbative results are 
generated through the use of a variational optimization pro- 
cedure known as the principle of minimal sensitivity (PMS) 
10]. The two dimensional Gross-Neveu model offers a per- 
fect testing ground for the LDE-PMS because, apart from 
sharing common features with QCD, it is exactly solvable in 
the large-A^ limit. The large-A^ result for the critical tempera- 



ture (at zero chemical potential) of the Gross-Neveu model is 
Tc ~ 0.567 m/r(0) where m/r (0) is the fermionic mass at T = 0. 
However, due to the appearance of kink-anti-kink configura- 
tions, the exact critical temperature for this model should be 
zero 17] . Because kink configurations are unsuppressed the 
system is segmented into regions of alternating signs of the 
order parameter, at low temperatures. Then, the net average 
value of the order parameter is zero. At leading order, the 
1 /N approximation misses this effect because the energy per 
kink goes to infinity as N ^ °° while the contribution from 
the kinks has the form e^^ . Our strategy will be twofold. 
First, we show that the LDE-PMS exactly reproduces, within 
the N limit, the "exact" large-A^ result. Next we show ex- 
plicitely that already at the first non trivial order the LDE takes 
into account finite A^ corrections which induce a lowering of 
Tc as predicted by Landau's theorem. Here, the calculations 
are performed for three cases which are: (a) T = and /j = 0, 
(h)T y^Q and /j — Q and (c) T = and /J 7^ 0. Our main results 
include analytical relations for the fermionic mass at T = 
and — 0, Tc (at fJ — 0) and /j^ (at T = 0) which include fi- 
nite A^ corrections. The case T ^0 and /J ^ 0, which allows 
for the determination of the tricritical points and phase dia- 
gram is more complex, due to the numerics. This situation is 
currently being treated by the present authors ]8]. In the next 
section we review the Gross-Neveu effective potential at finite 
temperature and chemical potential in the large-A^ approxima- 
tion. The LDE evaluations are presented in section III. The 
results are discussed in section IV while section V contains 
our conclusions. 



II. THE GROSS-NEVEU EFFECTIVE POTENTIAL AT 
FINITE TEMPERATURE AND CHEMICAL POTENTIAL IN 
THE LARGE-A' APPROXIMATION 
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The Gross-Neveu model is described by the Lagrangian 
density for a fermion field \\tk (k = I, . . . ,N) given by |5(] 



1 



N 
k=l 



(1) 



When nip = the theory is invariant under the discrete trans- 
formation 



(2) 



displaying a discrete chiral symmetry (CS). In addition, Eq. 
([U has a global SU{N) flavor symmetry. 

For the studies of the Gross-Neveu model in the large-A^ 
limit it is convenient to define the four-fermion interaction as 
^^A^ = X. Since g-^ vanishes like 1 /N, we then study the the- 
ory in the large-A^ limit with fixed A, [5J. As usual, it is useful 
to rewrite Eq. ([T]i expressing it in terms of an auxiliary (com- 
posite) field o, so that 



Wk{i ^)Wk~0\\tk\\ik- 



(3) 



As it is well known, using the 1 /N approximation, the large-A^ 
expression for the effective potential is 



--N 



2A. 



-iN 



cf-p 



In I 



(4) 



The above equation can be extended at finite temperature and 
chemical potential applying the usual associations and re- 
placements. E.g., momentum integrals of functions /(poiP) 
are replaced by 



/ig;/ta,P)-^.TX;/^/|,(»„-,>,),pl 



dp 



{2k) 



{2ii) 



where co,, = (2« + 1 )%T, « = 0, ± 1 , ±2, . . ., are the Matsubara 
frequencies for fermions fl^ . For the divergent, zero temper- 
ature contributions, we choose dimensional regularization in 
arbitrary dimensions 2co = 1 — e and carry the renormalization 
in the MS scheme, in which case the momentum integrals are 
written as 



+ T [ ln(l+exp{-[cOp(o,)+H/r}) 
Jp 

+ T [\n{\+t^^{-[<^ip{o,)-^i]/T}). (6) 
Jp 

After integrating and renormalizing the above equation one 
obtains 



A' 2A. 271 



where 



1 f M 



2T^h ia,b) 



(7) 



/ dx 
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Jo 





, (8) 



with a = (Tc/r and h—iijT. Taking the T = and /j = limit 
one may look for the effective potential minimum (d; ) which, 
when different from zero signals dynamical chiral symmetry 
breaking (CSB). This minimization produces ||5l@] 



m/7 (0) =0,. =Mexp(^-|) . 



(9) 




J (2n) ^ J„ 



(2k 



e/2 



471 / 



j2co. 
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i2(B 



where M is an arbitrary mass scale and Je — 0.5772 is the 
Euler-Mascheroni constant. The integrals are then evaluated 
by using standard methods. 

In this case, Eq. (HI can be written as 



yN 

^eff 



N 



21 



dp 



ln[(co„-/A')^ + «pK)] , (5) 



where (i)l{ac] 



The sum over the Matsubara's fre- 



quencies in Eq. ^ is also standard [10] and gives for the 
effective potential, in the large-A^ approximation, the result 



FIG. 1: A three dimensional graph showing the fermionic mass, nif, 
as a function of T and /j. One sees a second order phase transition at 
/J = while a first order transition occurs at J = 0. All quantities are 
in units of 10 x M while A. = jr. 



One may proceed by numerically investigating as a 
function of T and /j as shown in Figure [T| which shows a 
smooth phase (second order) transition at /j = 0. At this point, 
the exact value for the critical temperature (T^) at which chiral 
symmetry restoration (CSR) occurs can be evaluated analyti- 
cally producing ifTHl 



Tc = niF (0) — ~ 0.567 OTf (0) , 

7t 



(10) 



2 



Brazilian Journal of Physics, 37 ,258, (2007) 



25 
20: 
15- 
10 
5 



5 10 15 20 25 30 




5 10 15 20 25 30 

^ 



FIG. 2: Top views of figure 1. On the LHS we have a shaded fig- 
ure where the black region represents CSR. The contour lines of 
the figure on the RHS indicate an abrupt (first order transition) for 
small values of T. Both figures display a (tricritical) point where the 
smooth descent meets the abrupt one. All quantities are in units of 
10 X M while A- = Jt. 



where 11* = T| — (r| — 0^)5. The counterterm Lagrangian den- 
sity, Lct.5, has the same polynomial form as in the original 
theory while the coefficients are allowed to be 5 and r| depen- 
dent. Details about renormalization within the LDE can be 
found in Ref. |14^. 

From the Lagrangian density in the interpolated form, Eq. 
(fl4] i. we can immediately read the corresponding new Feyn- 
man rules in Minkowski space. Each Yukawa vertex carries 
a factor — /5 while the (free) o propagator is now —iX/ (NS). 
The LDE dressed fermion propagator is 



Sf{p) -~ 

where rj* = r] - (r| - 0^)5. 



IE 



(15) 



while, according to Landau's theorem, the exact result should 
be Tc — 0. By looking at Figure [T] one notices an abrupt (first 
order) transition when T — Q. The analytical value at which 
this transition occurs has also been evaluated, in the large-A^ 
limit, yielding [121 



wf(0) 

~7r 



(11) 



In the T — /J plane there is a (tricritical) point where the lines 
describing the first and second order transition meet. This can 
be seen more clearly by analyzing the top views of figure [T] 
Figure |2] shows these top views in a way which uses shades 
(LHS figure) and contour lines (RHS figure). The tricritical 
point (Pf^) values can be numerically determined producing 
Ptc = {TtcPtc) - [0.318mf (0),0.608mf (0)] M- 



III. THE LINEAR 5 EXPANSION AND FINITE A' 
CORRECTIONS TO THE EFFECTIVE POTENTIAL 

According to the usual LDE interpolation prescription 101 
the deformed original four fermion theory displaying CS 
reads 



N 
k=l 



(i ^) Va- + ri ( 1 - ?>)^ik\\ik + 5— i^^k^kf 



(12) 

So, that at 5 = we have a theory of free fermions. Now, the 
introduction of an auxiliary scalar field o can be achieved by 
adding the quadratic term. 



5A^ 



X 



2X y^ + N^^^"^^ 



(13) 



to L^{\\t,\\t). This leads to the interpolated model 



N 
k=l 



, (14) 






FIG. 3: LDE Feynman graphs contributing up to order-6. The black 
dot represents a 6t| insertion. The external dashed line represents Oc 
while the internal line is the a propagator. The last diagram brings 
the first finite A' correction to the effective potential. 



Finally, by summing up the contributions shown in figure [5] 
one obtains the complete LDE expression to order-5 



eff,5' 
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4n^N 1 



ln( — j ~l2ia,b) 

T 2 



In I — ) -/2(fl, 



-2T^Ii{a,l 



-Jl{a,b) kl6) 



where h{a,b) is defined by Eq. (O, with a — rj/r. Also, 



h{a,b) 



dx 



(b^-l 



and 



J2{a,l 



sinh(fe) 



dx- 



cosh(-\/?+fl^) + cosh(fe) 



(17) 



(18) 



Notice once more, from Eq. ( fT6l ), that our first order already 
takes into account finite corrections. Now, one must fix the 
two non original parameters, 5 and r|, which appear in Eq. 
JTSI l. Recalling that at 5 = 1 one retrieves the original Gross- 
Neveu Lagrangian allows us to choose the unity as the value 
for the dummy parameter 5. The infra red regulator r| can be 



3 



fixed by demanding V^g gi to be evaluated at tlie point where it 
is less sensitive to variations with respect to r|. This criterion, 
known as Principle of the Minimal Sensitivity (PMS) |@] can 
be written as 



eff,8l 



fi,5=l 



= 0. 



(19) 



In the next section the PMS will be used to generate the non- 



J 



perturbative optimized LDE results. 



IV. OPTIMIZED RESULTS 

From the PMS procedure we then obtain from Eq. ( fT6l ). at 
r| = fj, the general result 



'y{^,T,^l)+^—■yiy\J,^l) 



ri-o,. + r| 



2kN 



7CT^ d 
—m/T.,/T)-H^/T^,/T) 



= 0. 



Tl=fi 



(20) 



r 



where we have defined the function 



M 



^y^y\J,^i) = \n\^-yh{^|T,^i/T). (2i) 

Let us first consider the case N ^ °°. Then, Eq. ( |20l i gives two 
solutions where the first one is fj = Of which, when plugged 
in Eq. ( fTSI l. exactly reproduces the large-A^ effective potential, 
Eq. (|7]i. This result was shown to rigorously hold at any order 
in 5 provided that one stays within the large-A^ limit OlSll . The 
other possible solution, which depends only upon the scales 
M,T and jj, is considered unphysical ifisll . 



A. The case T = and = 
Taking Eq. ( |20] | at T = /j = one gets 



, M , 
In ( ^ I - 1 



rj-o, -r|:^^ln ( — 



'2%N 



M 



= 0. 



(22) 



As discussed previously, the first factor leads to the model 
independent result, fj = M /e, which we shall neglect. At the 
same time the second factor in (l22T i leads to a self-consistent 
gap equation for fj, given by 



The solution for fjgi obtained from Eq. (l23t is 



(23) 



iTlN Of. 



exp 



2kN 



r 2k.n 

Tjgi {Oc) = Mexp j + W 

(24) 

where W{x) is the Lambert W function, which satisfies 
W(x)exp[W(x)] =x. 

To analyze CS breaking we then replace rj by Eq. (|24] | in 
Eq. ( fTSI l, which is taken at T = and /j = 0. As usual, CS 



breaking appears when the effective potential displays minima 
at some particular value Of 7^ 0. Then, one has to solve 



^eff.5l(Coll-il5l 



5=l,o..=a, 



= 0. 



(25) 



Since mp — o^ , after some algebraic manipulation of Eq. 
and using the properties of the W{x) function, one finds 



nifiT ^ 0,iJ = 0)=MT (k,N) 1 



-1 



2N 



where we have defined the quantity J (A-,A^) as 



(26) 



f {X,N) = exp 



(27) 



X[1-1/(2A^)], 

Eq. ( |26] l is our result for the fermionic mass at first order in 



O-c 




X 



FIG. 4: The effective potential minimum,d(-, which corresponds to 
the fermionic mass, as a function of A. for A' = 1,3 and 10. The 
dashed line represents the large-A? result. Oc is given in units of M. 



5 which goes beyond the large-A^ result, Eq. (|9|l. Note that in 
the N °° limit, f (k,N — > 00) = exp(— Tt/A,). Therefore, Eq. 
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correctly reproduces, within the LDE non perturbative 
resummation, the large-A^ result, as already discussed. In Fig. 
|4] we compare the order-5 LDE-PMS results for Gc with the 
one provided by the large-A^ approximation. One can now 
obtain an analytical result for fj evaluated at Oc — Oc- Eqs. 
(l24l i and ^ yield 

r\5i{ac)=MfiX,N). (28) 

Fig. |5] shows that f[{Cc) is an increasing function of both 
and X kickly saturating for N >3. The same figure shows the 
results obtained numerically with the PMS. 




A 



FIG. 5: The LDE optimum mass (fj), evaluated at = Sc, as a func- 
tion of A. for A' = 1,3 and 10. The continuous lines were obtained 
from the analytical result, Eq. l|28), while the dots represent the re- 
sults of numerical optimization, fj is given in units of M. 



B. The case r / and ;U = 

Let us now investigate the case T y^Q and /j = 0. In prin- 
ciple, this could be done numerically by a direct application 
of the PMS the LDE effective potential, Eq. ( fT6l l. However, 
as we shall see, neat analytical results can be obtained if one 
uses the high temperature expansion by taking r[/T — a ^ 1 
and fj/T ^ b I. The validity of such action could be ques- 
tioned, at first, since r| is arbitrary. However, we have cross 
checked the PMS results obtained analytically using the high 
T expansion with the ones obtained numerically without using 
this approximation. This cross check shows a good agreement 
between both results. Expanding Eq. ^ in powers of a and 
b, the result is finite and given by lll6l] 



/l(fl< l,/7< 1) 



7t^ b^ a^ /7t\ a^ 

(29) 



and 



I,ia,b)^lnQ-j, + l^lb'+"-]+Oia\b')^ 



(30) 



where ^(3) ~ 1.202. If we then expand Eq. (fTST i at high tem- 
peratures, up to order rf-fT^, we obtain 



eff.5' 



21 6 



271: 



In 



l'Mey^\ 7^(3) Tl 



\ Tk J 4(27t)2r2 



,11(11 -a,) 



In 



/MeT£\ 7C(3) 11^ 
fMe^\ 7C(3) /M£^\ if' 



(31) 



Now, one sets 5=1 and apphes the PMS to Eq. {T\\ to obtain 
the optimum LDE mass 



fi(o,,r) = o, 1 



N(2%) 
7C(3) a? 



In 



2(27t)2 r2 



N(2n) 



]n 



Tn 



-2 



-1 



(32) 



The above result is plugged back into Eq. ( |3T1 ) which, for 
consistency, should be re expanded to the order ri^/T^. This 
generates a nice analytical result for the thermal fermionic 
mass 



Oc{T) = ± 



A?Vl47iC(3)A. 



2Nk + In 



-,3/2 



V Tk 



-2Nk+{2N- l)Xln 



Tn ) 



1/2 



(33) 



Figure|6]shows Oc{T) /M given by Eq. (|33]) as a function of 
T /M, again showing a continuous (second order) phase tran- 
sition for CS breaking/restoration. The numerical results illus- 




FIG. 6: The effective potential minimum, Oc, as a function of the 
temperature. Both quantities are in units of M and have been plotted 
for N = 3 and X = Tt. The dotted line corresponds to the large result 
predicting Tc = 0.208M while the continuous line, which represents 
the LDE result, predicts = 0. 170M. In both cases the transition is 
of the second kind. 



5 



trated by Fig. |6]show that the transition is of the second kind 
and an analytical equation for the critical temperature can be 
obtained by requiring that the minima vanish at Tc- From Eq. 
( |33] | one sees that (T = T^) = can lead to two possible 
solutions for Tc. 




X 



FIG. 7: The critical temperature (Tc ), in units of M, as a function of X 
for N = 1,3 and 10. The continuous lines represent the LDE results 
while the dotted line represents the large-A' result. 



The one coming from 



2Af7l + ln 



~T^) 



= 0. 



(34) 



can easily be seen as not been able to reproduce the known 
large-A^ result, when N ^ °°, Tc = Mexp(Y£ — k/X)/k. How- 
ever, the other possible solution coming from 



with fx{N) as given before, by Eq. (l27b . Therefore, Eq. 
also exactly reproduces the large-A^ result for N °°. The 
results given by this equation are plotted in Fig. [T] in terms 
of A- for different values of A^. The (non-perturbative) LDE 
results show that Tc is always smaller (for the realistic finite 
case) than the value predicted by the large-A^ approximation. 
According to Landau's theorem for phase transitions in one 
space dimensions, our LDE results, including the first 1 /N 
correction, seem to converge to the right direction. 



C. The case T = and a/ ^ 

One can now study the case T = 0,/j 7^ by taking the limit 
r ^ in the integrals I\ , I2 and J2 which appear in the LDE 
effective potential, Eq. ( fTSI l. In this limit, both functions are 
given by 



\xmTH,{a,b) = --Q{^i~^) 



n^ln(^+V^ 



limhia^b) = 6(u — r|)ln 



\lmTJ2{a,b) = ^gn{^l)Q{^l-^)^Jf~^ 



(37) 
(38) 
(39) 



Then, one has to analyze two situations. In the first, r| > /v, 
the optimized fj is given by 



= 0. 



(35) 



gives for the critical temperature, evaluated at first order in 5, 
the result 



Tc=M — exp<!— 

7t ^1 X[\-\/{2N)] 



= M—T{KN) , (36) 

7t 



Xi? 1 [M 

In — 

27tA^ T] \ T] 



ri-Oc 



11=11 



X^JM 
2%N \y\ 



(40) 



while for the second, r| < /j, r| is found from the solution of 



r| - a. 



A,r| 
27tA^' 



In 



M 



In 



M 



Xy\ 

2kN 



= 0. 



n=Ti 



(41) 



Note that the results given by Eqs. ( I37ll39] l vanish for fJ <r[- 
Fig. [S] shows fjc, obtained numerically, as a function of A- for 
different values of A^. Our result is contrasted with the ones 
furnished by the 1 /N approximation. The analytical expres- 
sions for fjgi (Oc), Eq. (|28] |. and Tc, Eq. ( l36b . suggest that an 
approximate solution for /Jc at first order in 5 is given by 



^lc{T^O)^^^(X,N) 



(42) 



It is interesting to note that both results, for Tc, Eq. (Q, and 
/Ac, Eq. ( |42] |. follow exactly the same trend as the correspond- 
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FIG. 8: The critical chemical potential /Jc in units of M, plotted as 
a function of X for N = 3 and T = 0. The dashed line represents 
the l/N result at leading order, the dot-dashed line represents the 
1 /N result at next to leading order and the continuous line is the first 
order LDE result. 



ing results obtained from the large-A^ expansion, Eqs. ( fTOl l 
and (fTTI) . respectively, which have a common scale given by 
the zero temperature and density fermion mass m/r(0). Here, 
the common scale is given by fj evaluated at Oc = Oc and 
r=/7 = 0,f|5i (a,)=Mj {X,N). 



dimensional Gross-Neveu model at finite temperature and 
chemical potential. Our results show that when one stays 
within the large-A^ limit the LDE correctly reproduces the 1/A^ 
approximation leading order results for the fermionic mass, 
Tc and /j^. However, as far as Tc is concerned the large-A^ 
predicts Tc ~ 0.567 mf(0) while Landau's theorem for phase 
transitions in one space dimensions predicts Tc — 0. Hav- 
ing this in mind we have considered the first finite correc- 
tion to the LDE effective potential. The whole calculation 
was performed with the easiness allowed by perturbation the- 
ory. Then, the effective potential was optimized in order to 
produce the desired non-perturbative results. This procedure 
has generated analytical relations for the relevant quantities 
(fermionic mass, Tc and /Jc) which explicitely display finite 
coiTections. The relation for Tc, for instance, predicts smaller 
values than the ones predicted by the large-A^ approximation 
which hints on the good convergence properties of the LDE in 
this case. The LDE convergence properties in critical temper- 
atures has received support by recent investigations concerned 
with the evaluation of the critical temperature for weakly in- 
teracting homogeneous Bose gases [17]. In order to produce 
the complete phase diagram, including the tricritical points, 
we are currently investigating the case r T^Oand^T^O |8|]. 
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